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THE WAVELET ELEMENT METHOD 
PART II: REALIZATION AND ADDITIONAL FEATURES IN 2D AND 3D * 

CLAUDIO CANUTOt, ANITA TABACCO*, AND KARSTEN URBANS 

Abstract. The Wavelet Element Method (WEM) provides a construction of multiresolution systems and 
biorthogonal wavelets on fairly general domains. These are split into subdomains that are mapped to a single 
reference hypercube. Tensor products of scaling functions and wavelets defined on the unit interval are used 
on the reference domain. By introducing appropriate matching conditions across the interelement boundaries, 
a globally continuous biorthogonal wavelet basis on the general domain is obtained. This construction does 
not uniquely define the basis functions but rather leaves some freedom for f ulfilling additional features. 

In this paper we detail the general construction principle of the WEM to the ID, 2D and 3D cases. 
We address additional features such as symmetry, vanishing moments and minimal support of the wavelet 

functions in each particular dimension. The construction is illustrated by using biorthogonal spline wavelets 
on the interval. 

Key words, wavelet element method, matching conditions. 

Subject classification. Applied and Numerical Mathematics 

1. Introduction. The construction of multiresolution systems and wavelets on general domains and 
manifolds in R n is a crucial issue for applying wavelet methods to the numerical solution of operator 
equations such as partial differential and integral equations. This problem has been recently addressed by 
many authors [15, 7, 13, 4, 14, 10, 16]. 

In [4], the Wavelet Element Method (WEM) was introduced borrowing ideas from analogous constructions 
in spectral methods. Tensor products of scaling functions and wavelets on the unit interval are mapped to 
the subdomains in which the original domain is split. By matching these functions across the interelement 
boundaries, globally continuous biorthogonal wavelet systems are obtained, which allow the characterization 
of certain function spaces and their duals. These spaces contain functions which are piecewise regular with 

respect to a Sobolev or Besov scale in each subdomain, and satisfy suitable matching conditions at the 
interfaces. 

The construction of [4] does not uniquely determine one particular wavelet basis. It rather indicates 
the algebraic conditions to be satisfied in order to obtain a globally matched basis, leaving the freedom of 
enforcing additional features for the basis. The purpose of the present paper is to exploit this freedom, 

-This research was partially supported by the National Aeronautics and Space Administration under NASA Contract No. 

n^l 704 *' WhUe thG firSt aUth ° r W8S ” residence at the Institute for Computer Applications in Science and Engineering 
(ICASE), NASA Langley Research Center, Hampton, VA 23681-2199. Other support come from the following funds- in 
Italy, MURST ex-40% Analisi Numenca, Progetto Strategic CNR Modelli e Metodi per la Matematica e PIngegneria, Progetto 
Coordinate) CNR Wavelets , analisi del segnale e simmetne and Progetto Vigoni ; in Germany, DAAD ViGONi-Project Multilevel - 
Zerlegungsverfahren fiir Partielle Differentialgleichungen and Volkswagen- Stiftung within the project Reduzierte Modellterung 
und Simulation von Vielstoffprozessen mit Multiskalen Verfahren. 

tDipartimento di Matematica, Politecnico di Torino, Corso Duca degli Abruzzi 24, 10129 Torino, Italy, and Istituto di 
Analisi Numerica del C.N.R., 27100 Pavia, Italy, e- mail: ccanutoepolito.it. 

t Dipartimento di Matematica, Politecnico di Torino, Corso Duca degli Abruzzi 24, 10129 Torino, Italy e-mail- 
tabaccoQpolito.it. ’ 

Slnstitut fur Geometric und Praktische Mathematik, RWTH Aachen, Templergraben 55, 52056 Aachen, Germany e-mail- 
urbanQigpm . rwth-aachen . de. 


1 



indicating how to construct wavelet bases that have properties such as, e.g., minimal local support, moment 
conditions or symmetry. 

While [4] deals with the general spatial dimension n, here we detail our construction in the univari- 
ate, bivariate and trivariate cases separately, showing additional properties in each particular dimension. 
Following the guidelines given in this paper, one can concretely build scaling functions and wavelets (as 
well as the related transforms) on fairly general domains starting from a given multiresolution analysis on 
the unit interval. As an example, we illustrate the results of our construction when a biorthogonal spline 
wavelet system on the interval is used. The filter coefficients for this example are given in the Appendix. An 
application to the numerical solution of 2D elliptic boundary value problems is presented in [5]. 

The outline of the paper is as follows. In Section 2 we recall the main properties of biorthogonal systems 
on cubes, which may fulfil homogeneous boundary conditions. Section 3 is devoted to the description of the 
multiresolution analyses on general domains. The construction of the corresponding wavelets in ID, 2D and 
3D is described in Sections 4, 5 and 6, respectively. 

2. Biorthogonal systems on reference domains. In this section, we use tensor products to con- 
struct scaling functions and biorthogonal wavelets on hypercubes starting from suitable multiresolution 
analyses on the unit interval (eventually fulfilling homogeneous boundary conditions). 

We will frequently use the following notation: by A < B we denote the fact that A can be bounded 
by a multiple constant times B, where the constant is independent of the various parameters A and B may 
depend on. Furthermore, A < B < A (with different constants, of course) will be abbreviated by A ~ B. 

2.1. General setting on the interval [0, 1). There are many examples of biorthogonal wavelets on 
the interval available in the literature, see [2, 9, 18, 6, 11, 17] for example. In this subsection we collect 
the main properties of those biorthogonal wavelet systems on the interval constructed in [11, 17]. We first 
describe the general approach and then the modifications for f ulfill i n g boundary conditions as introduced in 
[4, 13]. 

The starting point are two families of scaling functions 

3j := {t 3 ,k : k G Aj}, Zj := {&* : k G A,} C L 2 (0, 1), 

where A., denotes an appropriate set of indices and j > jo can be understood as the scale parameter (with 
some j 0 denoting the coarsest scale). For subsequent convenience, these functions will not be labeled by 
integers as usual, but rather by a set of real indices 

(2.1) Aj := {tj, 0 = T i , 1 <r ij 2 < - <T j)Ki = l- 

In other words, each basis function is associated with a node, or grid point, in the interval [0, 1]; the actual 
position of the internal nodes r J( 2» • • • . T. j . k , - 1 will be irrelevant in the sequel, except that it is required that 
Aj C Aj+i (see (2.3-h)). It will be also convenient to consider Ej as the column vector (£j,k)ke Aj, and 
analogously for other sets of functions. 

The families Ej , Ej can be constructed to be dual generator systems of a multiresolution analysis in 

L 2 ( 0 , 1 ) 

(2.2) Sj : = span Ej , Sj := spanE j? 
in the sense that the following conditions in (2.3) are fulfilled: 
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(2.3-a) The systems Ej and Ej are refinable, i.e., there exist matrices such that 

~j = Mj — j+i , s.j = Mj — j+i ■ 

This implies, in particular, that the induced spaces Sj, Sj are nested, i.e., Sj C S j+1 , Sj C S j+1 . 

(2.3-b) The functions have local support: diam (supp ~ diam (supp£j,jt) ~ 2~ J . 

(2.3-c) The systems are biorthogonal, i.e., fe,k, 4 *>')l 2 (o,i) = h,k>, for all k,k' £ Aj. 

(2. 3-d) The systems Ej, Ej are uniformly stable, i.e., 

T. ~ l|c||f 2 ( Aj) ~ y Ck£j,k 

keAj L*( 0 , 1 ) feeA, 1 , 2 ( 0 , 1 ) 

where c := (c k ) k€ A, - 

C 2 - 3 - 0 ) The functions are regular, i.e., £ jtk £ H"t{ 0 , 1 ), £ j<k £ H^( 0 , 1 ), for some 7,7 > 1, where H s ( 0 , 1 ), 
5 > 0, denotes the usual Sobolev space on the interval as defined , e.g in [1]. 

(2.3-f) The systems are exact of order L, L > 1 , respectively , i.e., polynomials up to the degree 
are reproduced exactly: P L -i(0,l) C S jy P L _ X ( 0,1) C S jf where P r (0, 1) denotes the set of the algebraic 
polynomials of degree r at most, restricted to [ 0 , 1 ]. 

(2.3-g) There exist biorthogonal complement spaces Tj and Tj such that 

Sj+i = Sj 0 Tj, Tj 1 Sj , §j+! = Sj © Tj, Tj 1 Sj. 

(2.3-h) The spaces Tj and Tj have bases 

~ {Vj^h : h € Vj}, Tj = {f}j,h : h 6 Vj}, 

(with Vj := A J+ i \ Aj = ■ . . , ^j,Mj }j 0 < Vj, i < •■■ < Vj,M j < 1) which are biorthogonal (in the 

sense of ( 2.3-c )) and uniformly stable (in the sense of (2. 3-d)/ These basis functions are called biorthogonal 
wavelets. 

(2.3-i) The systems Ej and Tj are boundary adapted, i.e., at each boundary point: 

(i) only one basis function in each system is not vanishing; precisely , 

( 2 * 3 ' i * 1 ) £j,k(0)^0 <=> k = 0 , £j,fc( 1)^0 A: = 1 , 

( 2 * 3 ' L2 ) Vj,h(°) ^ 0 <=> h = Vj,u VjJ i(l) 7 ^ 0 <*=► h = v jM .\ 

(ii) the nonvanishing scaling and wavelet functions take the same value; precisely, there exist constants 
Co and c\ independent of j such that 

(2.3-L3) 6,o(0) = in„. x (0) = Co2^ 2 , £m(1) = (1) = Cl 2^ 2 . 

The same holds for the dual systems Ej and Tj . 

(2.3-j) The system Ej is boundary symmetric, i.e., £j,o(0) = £j,i(l) =: Aj and also for Ej. □ 

In addition, thanks to Jackson and Bernstein type inequalities, these systems yield norm equivalences 
for a whole range in the Sobolev scale: 

00 2 00 

(2A) y c jo, k £jo, k + y y d M ~ y 2 2 * jo |c, 0 ,fci 2 + y 2 2sj \dj, h \ 2 , 

fe6A w j=joheVj X’ fee a jo j=j 0 fegVj 
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Fig. 2.1. Primal scaling functions. 

where s G ( - mm(L, 7 ),min(L, 7 )) and X s = H s ( 0, 1) if s > 0 or X s = {H~ B { 0, 1))' if s < 0. In particular, 
for s = 0we have a Riesz basis of L 2 (0, 1). 

The following concept will be important in the sequel. The system zz 3 is said to be reflection invariant, 
if Aj is invariant under the mapping i^l-i and 

(2.5) £ iik (l -a:) =£j,i-fc(a;), for all x e [0, 1] and k e Aj . 

A similar definition can be given for the system T j, as well as for the dual systems. If zij is reflection 
invariant, then Tj can be built to have the same property. This will be always implicitly assumed. 

Example . Throughout the paper, we shall illustrate our construction of matched scaling functions and 
wavelets starting from biorthogonal spline wavelets on the real line, as introduced in [ 8] . The corresponding 
multiresolutions on the interval are built as in [11, 12] with the choice of parameters L = 2 and L = 4, using 
SVD for the biorthogonalization. The particular implementation used to produce the pictures of the present 
paper is described in [3j. Figures 2.1 and 2.2 show the primal and dual scaling functions which are boundary 
adapted, whereas Figures 2.3 and 2-4 refer to primal and dual wavelets. These, and all the subsequent Figures 
of the paper, correspond to the level j = 4. “ 

2.2. Homogeneous boundary conditions on the interval. Boundary adapted generator and wa- 
velet systems can be easily modified to fulfill homogeneous Dirichlet boundary conditions. To this end, let 
us first introduce the following sets of the internal grid points: 

(2.6) A* nt := A j \ {0, 1}, Vf* := V; \ 

Let us collect in the vector 0 = (0o,Pi) € {0, l} 2 the information about where homogeneous boundary 
conditions are enforced, i.e., 0d = 1 means no boundary condition, whereas 0 d = 0 denotes boundary 
condition at the point d € {0, 1}. Correspondingly, let us set 

A* nt , if /? = (0,0), 

Aj \ {0}, if /3 = (0,1), 

Aj \ {1}, if /3 = (1,0), 

A j, if/? = (1,1). 

Let the generator systems be defined as 

=- 0 3 ■= {&,k ; k ^ A j }- : = (4* : k e A j }> 
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Fig. 2.3. Primal wavelets. 


and let us define the multiresolution analyses 


(2-8) Sj '.= span , S? := spanH^. 

Note that we have simply omitted the scaling functions which do not vanish at those end points of the 
interval where boundary conditions are enforced. 

The associated biorthogonal wavelet systems are the same as the previously defined ones, except 

that we possibly change the first and/ or the last wavelet depending on /?. More precisely, the wavelets can be 
chosen to vanish at each boundary point in which the corresponding component of (3 is zero. If the boundary 
condition is prescribed at 0, the first wavelets 7 ) jyl/j l and f]^ Ujl are replaced by 


(2.9) 





^2 £j, o)j 


respectively. The wavelets and rjP„ vanishing at 1 are defined similarly. Observe that the set of 

grid points Vj which labels the wavelets does not change, i.e., = Vj for all choices of f3. 
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Fig. 2.5. Primal and dual wavelet corresponding to the first wavelet grid point having homogeneous boundary conditions. 


The new systems E?, and s£, fulfill the conditions in (2.3) stated above, provided the index 
{3 is appended to all symbols. In addition, if the systems — j and are reflection invariant, see (2.5), the 
systems with boundary conditions can be built to be reflection invariant as well, in an obvious sense (i.e., 
the mapping x 1 - x induces a mapping of E^ into itself if /3 = (0, 0) or 0 = (1, 1), while it produces an 
exchange of with E^ 1 ’^ in the other cases). 

Example (continued). Figure 2.5 shows the modified wavelets defined in (2.9) for the B spline mul- 
tiresolution chosen to illustrate our construction. ^3 

2.3. Tensor products. Multivariate generators and wavelets can be easily built from univariate ones 
using tensor products. Hereafter, wc describe the construction in the domain H = (0, l) n , which will serve 
as a reference domain later on. 

Let the vector b = (J } 1 , . . . ,/? n ) (where (3 l 6 {0, l} 2 for 1 < l < n) contain the information on the 
boundary conditions to be enforced. Let us set, for all j > 

vfin) ~sf ®---®sf, 

and similarly for These spaces are trivially nested. In order to construct a basis for them, we define 

for x = (xi, . . . ,x n ) £ ft and k = (k \, . . . , k n ) € A j := Aj x - • • x Aj 

n 

~ ® ■ ■ • ® €>,*,)(*) = IIW**); 

1=1 

we set 

h ■■= {^,/t : * e A $}> 
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so that 


V/iQ) = span (*j), Vf(Cl) = span (Sj). 

Biorthogonal complement spaces Wj(fl), i.e., spaces satisfying 

Vri («) = v 3 b m © W}(b), w*((l) ± V/(n), 

are defined as follows. Set Vj := A J +1 \ Aj and for any h = (hi, . . . ,h n ) € V*, define the corresponding 
wavelet 

n 

1=1 

where 

i. . = /U' if e . 

Kfc. 

Then, set 

:= {^j,h : ^ e Y$}i and Wj(fi) := span^j. 

The definition of ^ and W^(fl) is similar. The wavelet systems and tyj are biorthogonal and form Riesz 
bases in L 2 (fl): the norm equivalences (2.4) extend to the multivariate case as well. 

Finally, considering the boundary values, we note that, given any l e {1, . . . , n} and d G {0, 1}, we have 

( < Pj ] k)\x l =d = 0 iff 7 ^ d or (ki = d and /3 l d = 0) 

and 

= 0 if f hi^i/d or (hi = v d and (5 l d = 0), 
with if d = 0, and z/ d = i /j yMj if d = 1. 

3. Multiresolution on general domains. In this section, we describe the construction of a mul- 
tiresolution analysis on our domain of interest ft C JR n . As already mentioned, fl is split into subdomains 
Qi, which are images of the reference element fl = (0, l) n under appropriate parametric mappings . The 
multiresolution analysis on fl is then obtained by transformations of properly matched systems on fl. 

We will first describe our technical assumptions on ft and the parametric mappings from the reference 
domain to the subdomains. In Subsection 3.2, we detail the construction of scaling functions and wavelets 
on a subdomain. The matching for the scaling functions will be described in Subsection 3.3. Most of 
the wavelets in ft simply arise by a mapping from the reference element without any matching; these 
functions will be described in Subsection 3.4. The more complex matching for the wavelets will be detailed 
in the subsequent sections for the ID, 2D and 3D cases, separately. Finally, in Subsection 3.5, certain 
characterization properties of Sobolev spaces are summarized. 

3.1. Domain decomposition and parametric mappings. Let us consider our domain of interest 
f} C JR n , with Lipschitz boundary dO^The boundary is subdivided in two relatively open parts (with 
respect to the Dirichlet part r D i r and the Neumann part r Neu , in such a way that 

dfl = f Dir U f Neui F Dir H T^ eu = 0. 

x Note that also certain domains that do not have Lipschitz boundary can be treated by this approach. Indeed, we have to 
assume that can be split into subdomains in such a way, that the assumptions on the parametric mappings are fulfilled. 
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We assume that there exist N open disjoint subdomains Sli C SI (i — 1, . . . , N) such that 

N 

n = \Jsi t 

*=i 

and such that, for some r > 7 (see (2.3-e)), there exist r-time continuously differentiable mappings F t : SI 
Cli (i = 1, . . . , N) satisfying 

(3.1) n t = Fi(fi), \JFi\ := det(JFi) > 0 in 6, 

where JFi denotes the Jacobian of F*; in the sequel, it will be useful to set Gi := Ff 1 . 

Let us first set some notation, starting with the reference domain SI. For 0<p<n-l,a p face of SI is 
a subset a C dSl defined by the choice of a set C a of different indices l \, . . . ,Z n - P G { 1 , . • . , n} and a set of 
integers d \, . . . ,d n _ p G {0,1} in the following way 


(3.2) c 7 — { (x \ , . . . , x n ) : x^ — d \ , . . . , x/„_ p — d n —p, and 0 ^ x; ^ 1 if l C , a } 


(thus, e.g., in 3 D, a 0 face is a vertex, a 1 -face is a side and a 2 -face is a usual face of the reference cube). 
The coordinates x t with l e C a will be termed the frozen coordinates of a, whereas the remaining coordinates 
will be termed the free coordinates of <r. 

Next, we set up the technical assumptions for the mappings F*. To formulate these, we need some 
notations. Let a and 0* be two p - faces of SI, and let H : a — ► a f be a bijective mapping. We shall say that 
H is order -preserving if it is a composition of elementary permutations (5, t) >— ► ( t , s) of the free coordinates 
of a. An order-preserving mapping is a particular case of an affine mapping (see [ 4 ], Lemma 4 . 1 ). 

The image of a p-face of SI under the mapping Fi will be termed a p-face of Sli\ if T * y := n dlV is 
nonempty for some i ^ i r , then we assume that is a p-face of both Sli and SI */ for some 0 < p < n — 1 . 
In addition, setting T iy = Fi(d) = F»'(a'), with two p-faces a and a f of St, we require that the bijection 
Hiy := Gi> o Fi : a — ► a* fulfills the following Hypothesis (3.3): 

a) H iy is affine; 

b) in addition, if the systems of scaling functions and wavelets on [ 0 , 1 ] are not reflection invariant (see 
( 2 . 5 )), then is order preserving. 

Finally, the decomposition is assumed to be conformal in the following sense: the intersection fit Pi Sl t ' 
for i i r is either empty or a p face, 0 < p < n — 1 ; moreover, for i = 1 , . . . , N we suppose that dSli O Toir 
and dSli D f Neu are (possibly empty) unions of p-faces of Sli. 

To summarize, we need the following assumptions: 

(a) The mappings Fi : SI — ► fi* are r-times continuously differentiable, with r > 7, and satisfy | JFi\ > 0 
in SI ; 

(b) The mappings F* fulfill Hypothesis 3.3 above; 

(c) The domain decomposition is conformal. 

3 . 2 . Multiresolution and wavelets on the subdomains. Let us now introduce multiresolution 
analyses on each 0 *, i = 1 , . . . , AT, by “mapping” appropriate multiresolution analyses on SI. To this end, let 
us define the vector b(Sli) ~ (/3 1 , . . . , /3 n ) G {0, l} 2n as follows 




0 , if F t ({xi = d}) C r Din 

1 , otherwise, 


/ = 1 , . . . ,n, d — 0, 1. 
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Moreover, let us introduce the one-to-one transformation v »— ► v := v o which maps functions defined in 
Cl into functions defined in Cl. Next, for all j > j 0) let us set 

VS(n<) := : v G v/ (n<) (ft)}. 

A basis in this space is obtained as follows. For any k G A^ n '\ set k 1,1 ' 1 := Fi(k); then, define the set of grid 
points 

AC} := {fc (i) : k G Aj (n,) } 

in Cli. The grid point A:, whose image under G t is k, is associated with the function in Vj(Qi) 

(t) * „ 

•Pj,k ~ ViJk ° 

i.e., (p^ k = tfj k . The set of all these functions will be denoted by 4>}. This set and the dual set <£} form 
biorthogonal bases of Vj (11,) and Vj ( O , ) . respectively, with respect to the inner product in L 2 (V„, ) 

(3.4) (u,v) := I u(x) v(x) | JGi(x)\ dx = / u(x)v(x)di :, 

which, due to the properties of the transformation of the domains, induces an equivalent L 2 type norm 

IMl!»(n,) ~ ( v , v )n t = ||«Hta ( n)' VveL 2 (fli). 

Coming to the detail spaces, a complement of Vj(Cli) in Vj+i(Hi) can be defined as 

WiW)i={v> ■ 

A basis ^ in this space is associated with the grid 

(3.5) U) := AC} +1 \ AC} = {h = F«(ft) : h G V?” 0 } 
through the relation 

(3.6) ■= Vh € H), h = Fi(h). 

The space and the similarly defined space form biorthogonal complements; the bases 

and are biorthogonal (with respect to {•, - Jn,). It is easily seen that the dual multiresolution analyses on 
Cli defined in this way inherit the properties of the multiresolution analyses on Cl. 

Finally, we introduce a concept that will be useful in the sequel. A point h € 7ij is termed internal to 
if h = with h = (hi, • * , h n ) such that each component hi belongs to U V* n * (see (2.6)). 

3.3. Multiresolution on the global domain. Now we describe the construction of dual multireso- 
lution analyses on Cl. Let us define, for all j >3o, 

(3-7) V,(n) := {w € C°(fi) : ^ e Vj^), i = l, . . . , N}- 

the dual spaces Vj(Q) are defined in a similar manner. In order to define a basis of V^ft), let us introduce 
the set 

N 

(3.8) K.j := (J AC} 
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containing all the grid points in ft. Each point of K,j can be associated to one single scale basis function of 
Vj(D), and conversely. To accomplish this, let us set 

I{k) := jt € {1, . . . , N) : k € 1%}, Vk e K,, 

as well as 

i fc (<) := Gi(k), Vi e I{k), Vk € Kj. 


Then, for any k G fCj let us define the function <pj t k as follows 


(3.9) 


<Pj,k in. 


UWI 1/2 if * € I(k), 

0, otherwise. 


This function belongs to Vj(D), since it is continuous across the interelement boundaries (see [4], Section 
4.2). Let us now set := {ipj,k ‘ k G /Cj}. The dual family := {<fj,h : k € tCj} is defined as in (3.9), 
simply by replacing each by <p ^ k . (For our B -spline example, matched scaling functions are displayed 
in Figures 4.1 and 5.2 in Sections 4 and 5, respectively.) Then, we have Vj(£2) = span4>j, Vj(Q) = span 
By defining the L 2 type inner product on Q 

N 

(3.10) ( u , v) n := ^(u, u) n . , 

t=l 

it is easy to obtain the biorthogonality relations from those in each fl*. 

3.4. Wavelets on the global domain. We now begin the construction of biorthogonal complement 
spaces Wj(Sl) and Wj(yi) ( j > jo) for Vj+i(Q) and Vj+i(ft) as well as the corresponding biorthogonal bases 
and Given the set of grid points 


N 

n, iCj+i \ iCj =: |Jw‘, 

i=l 


we shall associate to each h^Hj a function y: :J j L € and a function x pj t h € VV ; (T!). Then, we shall set 

: = {V»j -,h ■ heHj}, Vj := {ipj.h ■ h e Hj}. 

At first, let us observe that if h 6 H* is such that the associated local wavelet ipj (defined in (3.6)) 
vanishes identically on dfli \ dfl, then the function 


(3.11) 




rpj'lix), if x G Q it 
0, otherwise, 


will be the global wavelet associated to h. This situation occurs either when h is an internal point of fl* 
(recall the definition of internal point given at the end of the previous subsection), or when ail non internal 
coordinates of h correspond to a physical boundary (see Figure 3.1). 

The remaining wavelets will be obtained by matching suitable combinations of scaling functions and 
wavelets in contiguous domains. In each of the three coming sections, we shall detail the construction of the 
univariate, bivariate and tri variate matched wavelets, respectively. Precisely, 

- in ID, wavelets are matched across the interface between two contiguous subdomains (i.e., subinter- 
vals), see Section 4; 

in 2D, wavelets are matched around the common vertex of several subdomains (cross point), or 
across the common side between two subdomains, see Section 5; 
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Fig. 3.1. Scaling function grid points (circles) and wavelet grid points (crosses) in the subdomain Q t . The upper and right 
parts of the boundary of the subdomain belong to T. Wavelet grid points in the shaded area are associated to global wavelets 
constructed according to (3.11). 

in 3D, wavelets are matched around cross points, or around a common edge of several subdomains, 
or across the common face of two subdomains, see Section 6. 

Before going further on, let us make a general remark on vanishing moments. Wavelets on the reference 
domain satisfy the conditions 

I x r xp , ~ h (x ) dx — 0, Vh, V|r| < L — 1, 

Jn 

where x T = (x^ 1 , . . . ,x£ n ) and |r| = max t ri ; this follows from the fact that Vj(ft) contains the set ‘p£_ 1 (S7) 
of all polynomials of degree < L — 1 in each space variable. Unless a very special mapping is used, similar 
conditions in ft are not satisfied. However, they are replaced by analogous conditions, which still imply the 
compression property of wavelets. Indeed, Vj(ft) contains the subspace 

r L _ 1 (n) = {pec°(ti) : (Pinje *>£_!(«), Vi}, 

so that one has 

{p,*l>j,h)n = 0, VheHj, Vp g Vi_ x (Q). 

A dual condition holds for the dual wavelets Y>j,/i(x). 

3.5, Characterization of Sobolev spaces. At the end of our construction we shall obtain a system 
of biorthogonal wavelets on ft which allows the characterization of certain smoothness spaces. For instance, 
let us set 

(3.12) i/*(ft) := {v e H S (Q) : v = 0 on r Dir } 

for s 6 JN \ {0}, and let us extend the definition by interpolation for s £ iV, s > 0 (after setting H® (ft) = 
L 2 (Q)). Furthermore, we introduce another scale of Sobolev spaces, depending upon the partition V 
{fli : i = 1, . . . , AT} of ft; precisely, we set 

(3-13) := {v € Hl(Sl) : w,* € i = 1, . . . , AT} 

for s 6 JN \ {0}, equipped with the norm 

N 
i= 1 


li 




and we extend the definition using interpolation for s W, s > 0 (again we set = L 2 (£l)). The 

following Theorem summarizes the characterization features of our wavelet systems; they can be exploited 
in many different applications. 

Theorem 3 . 1 . (] 4 ], Theorem 5 . 6 ,) Assume that s E [ 0 , min(L, 7 )) . Then 

HZ(n-,V) = L € L 2 (n) : jr Y 2 2sJ l<v,^,A>n| 2 < oo 

( j~jo heHj 

In addition , if v E then 

00 

( 3 - 14 ) v= Y ( v i ( Pio,k)n<Pjo,k + Y { v ^hh)u^j,h, 

kefc j0 j=j 0 heHj 

the series being convergent in the norm of H£(Q; V), and 

00 

( 3 . 15 ) MUr-flw.) ~ Y 22 sJo l^^o,On | 2 + Y X 

k£Kj 0 j~jo /iCKj 


A dual statement holds if we exchange the roles 0 / V)(fl) and Vj(fl). 

Moreover , if s E ( — min(L, 7), 0) , the formulas (3.14) and (3.15) hold for all v E H£(Q;V) := 
(H^ (Q; V))' , provided the inner product (*, -)q is replaced by the duality pairing between the spaces H£(Q; V) 
and □ 

4. Univariate matched wavelets and other functions. In this section, we describe the construction 
of matched wavelets and other functions in the one-dimensional case. Since this material will be used in 
the subsequent construction of higher dimensional wavelets, we restrict ourselves to the natural reference 
situation of the interval I = (—1, 1) divided in the two subintervals /_ = (—1,0) and /+ = (0, 1) by the 
interface point C = 0. It is straighforward to reduce any other one dimensional matching to the present 
situation, by possibly introducing a suitable parametric mapping. 

The scaling function associated to the interface point C — 0, is defined by 


(4.1) 


<Pj, o{x) 


1 ( $j,i(£ + 1 ), xel-, 

%/2 Uj,o(x), x€/+. 


Example (continued). For our B spline example , the function (4.1) and its dual are displayed in 
Figure J^.l. □ 

4.1. Wavelets. Let us consider the local basis functions on each subdomain 

fj,i(x + l), e = 0, 


i’e (x) := 

4 >, 


e—1, 
e = 1 » 


xel~, 


and let us set (see [4], (5.6)) 


V#-i( J ±) span{^ : e E { 0 , 1 }}. 
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J\ 

1 ' ’ 1 1 ' ' ' ' 

\K : 

i_JL 

!■ ■ ■ - 1 ■ ■ ■ - ^ 


Fig. 4.1. Matched primal and dual scaling functions at the cross point. 


Any function v ± € V^ fl (/±) can be written as v± = 5Zee{o,i} a t ^ us denote by a ± := (aj , af ) the 
vectors of the local degrees of freedom. We want to build the local space 

(/):= {veC°(I) * V \J± ^ Vj + i(I±)} 

by matching functions in V® +1 (I„) and and we want to find a basis for the subspace Wj(I) := 

{v 6 V/ +1 (J) : (Uj^o) l 3 (/) = 0}. The matching condition between two functions v ± 6 V^ +l {I±) reads: 

a o "I" a i (1) = a o Cj,o(0) + ttj 1 " ^>^^(0), 

which, in view of the boundary values of the univariate scaling functions and wavelets (see (2.3-i.3)), is 
equivalent to 

(4.2) a Q -{-a 1 = aj + a^. 

Next, we enforce the additional condition of orthogonality to <fj 0 which is expressed as 

°= (<*o &,i(- + !) + (• + 1 >»^j,cr l[ _ 1 .oi) x>a( _ 1>0) 

( 4 - 3 ) + ( a o &,<>(•) + «i'^> j ,i(-).^,c ll o, 1J ) J . 2(0 1} 

= 7/=(«o +of), 

where the last equahty is a consequence of the biorthogonality on the interval. 

Using the matrix-vector notation of [4] (see formula (5.24) therein), condition (4.2) and condition (4.3) 
multiplied by \J2 read 


T>ol = 0 , 

where 

a = (a“,a + ) t = (qq ,aq ,aj,a^)*. 

It is easily seen that 'DT) t = diag (4, 2) is positive definite, so that T> has full rank 2. This implies that 
dim Wj(I) = 2. It is directly seen that 

(4.4) Ker T> = span {(0, 1, 0, 1)*, (1, -1, -1, l) f }. 
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For the dual system, we have the same condition. So, it remains to find 2 particular choices of a G KerX> 
and d £ Ker T>, i.e., 


a 1 := a, A ( 0 , 1, 0 , 1 )* + a ;, 2 ( 1 , - 1 ,- 1 , 1 )*, 
d' := 5j t i (0, 1,0, 1)* + a ( , 2 (1, -1, -1, l) f , 


that will define primal and dual wavelets / = 1, 2, as 


1) + (a M - (HI), x £ 

-aj,2£j,o(£) + Ou + £ £ /+. 


The coefficients have to be chosen in order to obtain biorthogonal functions. Using the biorthogonality on 
the interval, it is readily seen that 


Sl,m = 


2ai,id m j +4 aj t2 d m , 2 . 


This can be rephrased by the matrix equation 

(4.6) Id = AX A 1 , 


where 



Since we have 4 equations for 8 unknowns, one can, in principle, choose 4 coefficients and the remaining 4 
are then determined by (4.6). 

It is convenient to relabel the wavelets t/r , / = 1,2, obtained by any particular choice of the 

coefficients as 


(4.7) 


M>i := ty, ■= 


(and similarly for the dual wavelets) so that they are associated in a natural way to the grid points h± £ Hj 
located around C and defined as /i + := F+(i/j t i), h- := F_(— 1 + 

The above mentioned freedom in the construction can be used to f ulfill additional features, which will 
now be described. 


4.1.1. Additional features. Depending on the particular application one has in mind, one might need 
the basis functions to have some additional features such as, for example, zero values at the cross point, 
(skew-) symmetry and reflection invariance. We will now address these issues. 

Zero values at the interface. The scaling and wavelet systems on the interval are supposed to be boundary 
adapted and boundary symmetric. Then, (4.5) implies 

i> l j( 0) = A j ai } /, t/^(0) = A j di,/, for l — 1 or l = 2, 

where A j is defined in (2.3-j). Choosing these coefficients to be zero implies zero value of the corresponding 
wavelet function at the interface. Note that it is not possible to enforce the condition for l — 1 and l = 2 
simultaneously, as this would contradict (4.6). 
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Symmetry and skew symmetry. Let us assume that Ej, Ej and T^, T j are reflection invariant. This 
implies in particular 




Then, we get, for x e I, 


^(ai,a 2 ) _ f _a 2 ^j,l(l + *) + ( a l + a 2)Vj,i + *)i X e 

3 \a 2 ^,o(x) + (ai -a 2 )% I/i l (x), x e I+, 




(*)> 


where the notation ^// ai ’“ 2 > abbreviates that the function is associated to the vector (ai,o 2 ). This implies 


ipj-x) = 


_/-0}(x), if oi,i = 0, 


V’j(i), if ai,2 = 0, 


for / = 1 or l = 2 , 


so that we can choose one wavelet function to be either symmetric or skew-symmetric. The choice A = 

( 0 1 \ 

[ i 1 leads to one symmetric and one skew -symmetric wavelet. The dual wavelets are then defined by 
^ .4 / 0 1/4 \ 

the matrix A = ^ ^ J and have the same properties. 


Reflection invariance. Under the same assumptions as before, the choice A 


-C‘.) 


leads to 


primal wavelets that reflect into each other under the mapping x — x. The dual wavelets are then defined 

1/4 1/8 


by the matrix A 


-( 


, and have the same property. 


1/4 -1/8 

Vanishing moments. Set, for a function / and a domain D, 


M r (/;£>)= f x r f{x)dx. 
Jd 


Using the fact that the scaling functions in [0, l] reproduce the constants, and exploiting the biorthogonality 
and the boundary adaptation property, it follows that 

(UO,o)l2(/+) = A” 1 = + 1))l2(/_), 

where again AJ 1 is defined in (2.3-j); by similar arguments, we get 

i 0,o)l 2 (/+) = 0 = (x y£j,i(‘ ~h 1 ))l 2 (I- ) l^r<L — 1. 

It is easily seen that these relations imply that the monomials x r (0 < r < L — 1 ) restricted to the whole 
interval I belong to Vj(I). Since the matched wavelets are orthogonal to Vj(I), we conclude that 

M r (^j;/) = 0, 0 < r < L — 1. 

Dual relations hold for $ ! • . We conclude that the wavelet functions that arise by the matching procedure 
automatically have the same order of vanishing moments on the whole interval I as the original functions 
on the interval [ 0 , 1 ]. 
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Fig. 4.2. Matched primal and dual wavelets at the cross point. 


Support located in only one subdomain. It is readily seen that (4.5) implies that the arising wavelets are 
both located in both subdomains. It is not possible to localize them on only one side. However, since 

diam (suppr^j^) ~ diam (supp^/t) ~ 2~ 3 , 

one still has diam (suppi/^/i) ~ 2~ J for h = h±; the same property holds for the dual functions. 

Example (continued). For our B spline example, matched wavelets defined in (4.5) with the choice of 
matrices A and A which guarantee reflection invariance (see above), are shown in Figure J^.2. □ 

4.2. Another basis of matched functions. Now we aim at defining a basis of the local space V^+i (/). 
Compared to Subsection 4.1, the orthogonality condition (4.3) is missing, so that we obtain the following 
matching conditions in matrix vector form 

Va = 0, with V:= (1,1, -1,-1), a = (a _ ,« + ) t = (a„ ,af 


It is obvious that Ker C = span {(1, —1, 0, 0)*, (0, 0, 1,-1)*, (0, 1,0, l) 4 } and we have to find 3 particular linear 
combinations, i.e., 

a‘ = 01,! (1,-1, 0,0)* + 01,3 (0,0, 1,-1)* + ai,3 (0,1, 0,1)*, l = l23 

o' = aj,x (1,-1, 0,0) t + Oi, 2 (0,0,1, -l) f + ai, 3 (0,1,0, l) 1 , 

which give rise to the three basis functions 


(4.8) 


f a l, 1 ( a /, 3 (*^)» ^ ^ 1—1 

l a /,2 £j,o(£) + (^1,3 “ a i,2) (£)» X€l+ y 


In this case, the biorthogonality gives the conditions 


= 2(a^id mj i -f a/,2d m ,2 + cn^dm^) — aj,ia m ,3 ~ a/, 3 a m ,i — a/, 2d m , 3 — a*, 36m, 2, 
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which can be rewritten as Id — B Y B , where 


/ ai,i 

a l,2 

a l,3 ^ 



dl,l 

di,2 

di >3 

\ 

( 2 

0 -1\ 

B := 1 a2,i 

02,2 

«2,3 

,B:= 


02,1 

«2,2 

02,3 


, Y := 0 

2 “M 

V a 3,l 

0 3 ,2 

a 3,3 / 


\ 

d 3 ,i 

d 3 ,2 

03,3 

) 

V-i 

-1 2 ) 


After choosing one particular solution of this algebraic system, we relabel the functions as 
( 4 . 9 ) S- := tf], $:=%, 

the dual functions are defined similarly. 

4.2.1. Additional features. Additional features can be required to the functions i 9 *-, l = 1,2, 3 just 
introduced. 

Vanishing moments. In this case we have 

Af r (#*•;/) = (a/,i + a/,2) M r (£j^;I), 

so that one preserves the order of vanishing moments if and only if a/,i = —a/, 2- 

Zero values at the cross point . Similarly to the results of Subsection 4 . 1 , we obtain i9*(0) = a /,3 and 
^(0) = d/, 3 . 

Symmetry of the arising functions. Using the same arguments as above, we obtain 

i' a2 ’ Q3) (_i) = # a2 ' a 1 ’ a 3 )(x), 


so that we have 


tf'(-x) 


-#*(£), if a/,i = -a/,2 and a/, 3 = 0, 

ti l Ax), ifoj,i=aj, 2 . 


Support located in only one subdomain. Obviously, one has 

f if ay 2 

( 4 . 10 ) suppt?* C i 

J l /+, if a/,i 


if a/,2 = a/,3 = 0, 
a/,3 = 0, 


so that it is possible to construct 3 functions, such that only one of them is localized in both subdomains. 

Example (continued). We give one particular example of three functions, one located in one m 
J. and one in both subintervals. For the latter one, also the vanishing moment property is preserved . Let 




( i 

0 

o > 

( 4 . 11 ) 

B := 

0 

0 

1 



l 0 

1 

0 / 


It is readily seen that B 1 = B l , hence we obtain the coefficients for the dual functions by 



l 3/4 

1/4 

1/2 

( 4 . 12 ) 

B = BY~ l = 1/2 

1/2 

1 


V 1/4 

3/4 

1/2 


These particular functions are displayed in Figure 1^.3 . 
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Fig. 4.3. Primal and dual bcLsis functions defined by (4.11) and (4.12). 


5. Bivariate matched wavelets. We shall now construct matched two dimensional wavelets, by 
firstly considering an interior cross point, next a boundary cross point and finally the common side of 
two subdomains. In each case, we indicate how wavelets can be defined, which have the most localized 
support. 

5.1. Matched wavelets around an interior cross point. We describe the construction of wavelets 
associated to grid points h e Hj which are close to a cross point C, at which Nc subdomains meet. We 
assume that these subdomains are (re-)labeled by fti, . . . ,Otv c , in a counterclockwise order. Moreover, we 
set r <(<+1 := dQ{ . We start with the situation in which C is interior to f2. In this case, it is convenient 

to set j . — fi ] . 

At first, we deal with a particular choice of the mappings to the reference domain. Next, we shall show 
that all other possibilities can be easily reduced to this choice, which - therefore - can be thought of as a 
reference situation . 

So, suppose that, for all i e {1,..., N c }, one has C = F;( 0,0) and r M+1 = Fi(d 0 i) = F i+ i(<7 10 ), where 
0oi := {(0,x 2 ) : 0 < £ 2 < 1} and <7i 0 := {(xi,0) : 0 < x\ < 1}. The grid points surrounding C, to which we 
will associate the matched wavelets, are the 2 Nc points hc,i defined as follows: 

(^•1) hc, 2 i-\ ~ hc, 2 i = F*(0, 1 < i < Nc- 

Note that each evenly numbered point belongs to a side meeting at C, whereas each oddly numbered point 
is internal to a subdomain meeting at C (see Figure 5.1). 

Dropping the index j, let us set 

'Poo (x) = 'Poo(x) = £j,o(xi)tj,o(£2), 

1 Poi( x ) = 'Poi(x) = 
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Fig. 5.1. Grid points hc,i, i = 1, . . . , 10, around a cross point C common to 5 subdomains . 


V’lo (x) = 1,Mx) = (Xl) fj,o(x 2 ), 

V'n (x) = ^n(x) = l (ii) tl (x 2 ), 

(see [4], formula (5.5)), as well as := span{^ !) : e G £ 2 }. A function is written as 

« (,) = ^2 q£ } V4 l) ; 

e€F 2 


we shall introduce the column vector a M := (a^) eG£ 2 . In order to characterize the local space 
: = £ C°(fl) : V|n 4 6 if i e {1, . . . ,N C }, v\n t = 0 elsewhere}, 


we proceed as in [4] (see Section 5.2), i.e., we enforce the continuity among the i/ 1 ) by considering the point 
C firstly, and the sides r^+i, secondly. Recalling (2.34.3), the continuity at C yields the set of linearly 
independent conditions 


(5.2) 


><0 

*oo 


+ a « 

-h a 01 


+ Q (l) 
+ a io 


4- a 


.(0 - 

a oo 


+ a ( 


(•+i) 

01 


< a io < a ll ’ 


1 < i < N c - 1, 


or 


cq ■ = Co • a (t+1 \ 1 < i < Nc - 1, 


with Co = (1, 1, 1, 1). Denoting by a := (QW)i=i iiM| jv c € ]R 4Nc the column vector of all degrees of freedom, 
these conditions can be written in matrix vector form as 


( Co -Co ^ 

Co a = 0, where Co — 

\ Co -Co / 

Let us now enforce continuity along the sides r\ ii+1 . To this end, observe that 


6 ]R( N c- l ) x4N c ' 


(5.3) 


v 

V 


|r M+1 ( x ) 

(*+i) 

|r<,, 


,(x) 


A J f( Q oo + “io)Cj,o(x 2 ) + (aioi +an)’0j, l , j ,,(x 2 )), 

A j(( a oo + J + a oi +1) )Cj,o(xi) + («i* 0 +1) + a ii +1 *)%i'j,i(xi)), 
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with A j defined in (2.3-j). Because of the linear independence of the univariate functions, the matching is 
equivalent to 


(5.4) 


Q W 

a oo 

(i) 




- JA+ 1 ) -L JJ+ 1 ) 


Of, 


01 


- 


a 


x oo 

(i+l) 

10 


-f a 


01 
(<+ 1) 
li 


Since we have already enforced the continuity at C G rj+i (see (5.2)), it is enough to require that a 
particular linear combination of the latter equations holds; precisely, we enforce 


JO JO 


JO 


(0 


*00 


X 01 


+ *10 -<*\{ =«, 


*00 


+ «l 


,(*+1) 

01 


— a 


(i+i) u+i) 


10 


— a 


li 


(for details, see [4], Proposition 5.2). Introducing the vectors d = (1,— 1,1,— 1) and d f = (1,1,— 1,-1), 
these conditions can be rephrased as 


( d -d' 


(5.5) 


C i a = 0, where C\ 




V -c" 


C' ) 


We are interested in finding a basis for the subspace (fl) := {v G : (v, (pj,c)n = 0}- Recalling 

that, by (3.9), 


ftj,C |f2i i x ) 

we obtain the condition 


x = Gi(x), i = 1, . . • , N C) 


E 


! 00 


i=l 


= 0. 


So, introducing the vector b (1,0, 0,0), all conditions enforced so far can be summarized in 


X>a = 0, where 



^ jp2Nc x4 Nc 


8 :={b,...,b)e JR lxiNc . 


It can easily be seen that 



g ]fd^ c x * 


the evident symmetric positive definite character of 7>T> 1 means that T> has full rank. This implies the 
existence of exactly 2Nc linearly independent functions in i.e., dim Wj^fi) = 2Nc- The parallel 
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construction for W? (Q) leads to the condition Dot = 0 (note that the matrix D is the same as for the 
primal system). Finally, we enforce biorthogonality between the primal and dual basis functions determined 
in this way. To this end, we choose a basis {«i, . . . , k- 2AT c } in Ker^D, and we look for linear combinations of 
these vectors 


2 N c 2 N c 

Ot ^ ^ n i m K'm ? Ot — ^ ^ &l,m t 1 < l ^ 2./V Q , 

m=l m=l 

such that, if a* is decomposed as the corresponding wavelets 

f V] a£ (,) V>e(i), if X € S7i, 1 < i < AT C , 

[ 0, elsewhere, 

and the dual ones wj y c defined in a similar manner, form a biorthogonal system. Then, setting /C := 

(^m)m=l,.,,,2iVci • (^i, m )^ m=l,. . .,2Nc > ^ • — (n^rn)i,7n=l,...,2Nc f 3-Ud *4. ! = (Oi )/=l,...,2A^cr = J&A j «A 1 = 

(o: i )/=i 1 ... ) 2Nc = JC-A , and exploiting the biorthogonality property in each subdomain, we can express the 
biorthogonality condition in the form A. = Id, i.e., 

(5.6) AK'KA 1 = Id. 

Since K, obviously has full rank, JC l K is regular, so this matrix equation has a solution. In particular, one 
can choose A = Id and consequently A = {K l K)~ l . In Subsection 5.1.1, we shall exhibit a specific basis in 
KerZ>, which allows to obtain primal wavelets with minimal support around C. 

Once the biorthogonal wavelets ip l j C and ^ c (1 < l < 2 Nc) have been determined, they can be 
associated to the 2 Nc grid points hc,i surrounding C, defined in (5.1). 

Reduction to the reference situation. Let us now show that we can reduce any interior cross point 
situation to the one described above. To this end, let us consider any subdomain Q. j having C as a vertex. 
Then, we have the following 4 cases: 
a) C = Fi(0, 0), b) C = fi(l,l), 

c)C=F i (0,l), d)C = F t (l,0). 

Recalling the assumption det ( JFi) > 0 in (3.1), it follows that in cases a) and b) the indices of the frozen 
coordinates of Ti-i^ and are given by _ i . = {2}, £p. = {1}, whereas in cases c) and d) one 

has Cp ,_ 4 . = {1}, Cp i = {2}. It is straightforward to see that the matching conditions along and 

r M +i in cases a) and b) yield the same vectors d and d f defined above, whereas the roles of these vectors 
are interchanged in the remaining cases. 

If = (aoo, ttoi^ a io, a u) 4 » us denote by := («oo^ a io^ a on a n ^he modified vector and 
by D the matrix obtained by modifying D according to cases c) and d). Since it is readily seen that 
Co * (d) 1 — Co * (d f ) t — 0,6- (c')* = 6 ■ (c' / ) < = 0, d * acM = d ( * as well as d ' • = d * a^, we obtain, 

as desired, Dot = D a. 

This procedure can be applied to all subdomains meeting at C, and so we are back to the reference 
situation. 

The biorthogonalization is performed following the same guidelines described above; obviously, the 
definition of the wavelets and the associated grid points has to be adapted to the specific orientation of the 
mappings to the reference domain. 
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5.1.1. Wavelets with minimal support. The perhaps most important feature, in view of numerical 
applications, is a minimal support of the wavelets, since this implies minimal length of the corresponding 
filters. 

Let us first consider under which circumstances it is possible to construct wavelets that are supported 
in only one subdomain fl*. By the matching at the cross point and the sides, we obtain the three conditions 

(5.7) Co • = c' * = c" * = 0. 


It is readily seen that the vectors (a, —a, —a, a), for some a G JR, are the only solutions of (5.7). Thus, any 
function defined as 


(5.8) 


f aTpoo(x) - aipoi(x) - i0 (x) + aV>n(£), x G f 

1 0 , elsewhere, 


belongs to and is supported only in ft*. However, such a function cannot be a wavelet because, 

by imposing the orthogonality to one gets the extra condition b • = a^o = a = 0 . This implies 

that it is not possible to have wavelets supported in only one subdomain. On the contrary, we are going to 
show that it is possible to construct 2Nq — 1 out of 2Nc wavelets to be supported in only two contiguous 
subdomain. To this end, let us deal with subdomains and fij+ j. Considering the block structure of the 
matrix the corresponding vectors of coefficients a^, a^ +1 ) have to fulfill the local equation 


Co —Co 


(i) 



OC ’ 

c' -c" 


a (*+l) 

b b 


L - 1 


It is easily seen that the two linearly independent vectors 


(« M0 , q 1,( * +1) ) := (0, 0, -1, 1, 0, -1, 0, 1)‘, 


(a 2 ’( i ),a 2 ’(* +1 )) := (-1, 1,0,0, 1, -2, -1,2)* 


solve this equation. Setting all the remaining coefficients in a to zero leads to the two seeked wavelet 
functions. 

We observe that one out of these 2 Nc functions is linearly dependent on all the others: since all these 
functions vanish at C, there are at most ( 2 Nc — 1) linearly independent functions among them. On the 
other hand, we now show that the following 2 Nc functions are linearly independent: 




Y Q e’ W 'Pe(i), 

x G Hi, 

(5.9) 

II 

0 

e£E 2 

Y a e’ ( ‘ +1) V’e(x), 

x G fJ-j+i 5 



0 , 

elsewhere, 

for / = 2 i 

— 1, i € {1,.. .,Nc}, 

Y «e' (,) V'e(x), 

X G 

(5.10) 

A 

II 

e£E 2 

Y a 2 /' +1) rp e (x), 

x G 5 



e€ E 2 
0 , 

elsewhere, 

for l = 2 i, 

i € { 1 , . . . ,Nc ~ 1}, and 



(5.11) 


f '0ii(*^)> x G fii, 1 < i Y Nci 

1 0 , elsewhere. 
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Indeed, if V 7 ™ c = ^ en V 7 (^) = 0 implies A 2 n c ~ 0 and i/>\r Nctl = 0 implies A 2 at c _i = 0. 

Then 'ift\n 1 = 0 implies Ai ^j,c + ^2 which in turns implies Ai = A 2 = 0. Considering in sequence 

0 2 , . . . , fijv c -i shows that all the coefficients must vanish. 

Note that the primal basis we have just exhibited corresponds to choosing the following basis in Ker 'D 


Km 

:= (0,.. 



if m = 2i - 1, i e {1,... , N c }, 

Km 

“ (0... 

.,0, (a 2, W)*, 0, . . 

..,0)*, 

if m = 2i, i € {1,...,N C - 1}, 

«2 N c 

:= (<*,•■ 



with d = (0,0,0, 1), 


and to choosing A = Id in (5.6). Then the dual matrix A is given by A = (X^JC) -1 , where 


K l K = 


Example (continued). We consider the situation in which four subdomains fh, . . . , meet at C (for 
simplicity , we assume linear parametric mappings from each subdomain to the reference domain). For our 
B-spline example , we show in Figure 5.2 the matched scaling function and the three different types of wavelets 
associated to grid points around C . □ 

5.1.2. Tensor products of matched univariate functions. A common situation for an internal 
cross point is the case Nc = 4, i.e., four subdomains meeting at C . In such a geometry, it is easy to 
construct a basis for Wj' (f2) by properly tensorising the univariate matched functions defined in Section 4. 

First of all, let us note that, by possibly introducing appropriate parametric mappings, we can reduce 
ourselves to the situation in which each subdomain is the image of one of the subdomains x I± (we use 
here the notation set at the beginning of Section 4), and C is the image of C = (0, 0). Then, let us consider 
the set of univariate functions given by the scaling function (pj ) 0 defined in (4.1), the wavelets xjjj defined in 
(4.5) (4.7) and the functions defined in (4.8)-(4.9). A basis in Wf(Q) is obtained by taking the image 
of 8 linearly independent tensor products of such functions satisfying the condition of orthogonality to the 
dual scaling function (pj 0 = (fj Q 0 (pj Q . For example, a possible choice is 

i>j ® » <Pj , 0 ® V 7 /, 0 

® 1 < pj , 0 ® V 7 / 0 

(obviously, these functions are extended by zero outside the union of the four subdomains) whose association 
to the 8 wavelet grid points around C is self evident. Note that this construction does not necessarily require 
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Fig. 5.2. Matched scaling function (1st raw) and primal wavelets In the first row the wavelet ip^ c is displayed. 

The second row shows the functions ipj c and ifrj c that are supported in Sli and The remaining 5 wavelets are rotations 
of these two functions. (Note that only a portion of each subdomain around C is shown.) 


the functions to have minimal support in the sense of (4.10), although efficiency will be enhanced by this 
feature. 
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Fig. 5.3. Grid points h c ,i (labeled only by i) around a boundary cross point in the 3 different cases (from left to right: 
pure Neumann, mixed and pure Dtnchlet case). 


5.2. Matched wavelets around a boundary cross point. Let us now consider the situation in 
which C E d$l is common to Nq subdomains Jlj, . . . , ordered counterclockwise. We assume that 
{^Ncn resp.) has a side, termed Tj (Fjv CJ resp.), which contains C and lies on d£l. Then, the following cases 
may occur: 

a ) Ti, t Nc e r Neu 

b) Ti E Toir, Tjg c E r^eu c) Ti E r ]\j eu , r Nc E r Dir 

d) Ti, T Nc e Toir- 

Since all arguments concerning biorthogonalization carry over from the interior cross point case, in the sequel 
we will only detail the matching and orthogonality conditions for each case separately. 

Pure Neumann case. Let us start by considering case a). The matching at the cross point C only differs 
from the interior cross point case by the absence of a matching condition between and Ct Nc . However, 
in the interior case, this condition turned out to be linearly dependent on the other ones, hence, there was 
no need to explicitly enforce it. This implies that the matrices describing the matching at C are the same: 

C Neu /» 

0 — ^ 0 * 

The matching conditions along the sides are the same as well, with the only difference that now the last 
row in (5.5) is missing: 


£»Neu 


d -c" 

\ c' — c" 


€ J?( w c7-l)x4JVc 


Since C £ I^Neui observe that there exists a dual scaling function associated to C . Consequently, we 
have to enforce orthogonality to this function, so we end up with the set of conditions 


Z> Neu a = 0, with X> Neu 


Co 

£*Neu 

B 




As above, it is easily seen that X> Neu has full rank; this implies that dim Wf(Q) = 2N C + 1. This is precisely 
the number of grid points surrounding C to which these wavelets can be associated; in the reference situation, 
they are the points (5.1) and the point h c , o := F\ {v jA , 0) (sec Figure 5.3, left). 

Mixed Neumann/ Dirichlet case. Obviously, the mixed Neumann/Dirichlet cases b) and c) can be viewed 
as symmetric ones, so we will only detail case b) here. Since C € lD, r , all functions in V)^j(Q) have to 
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vanish at C. This means that we have to add one more condition to those posed at C, i.e., we obtain the 
matrix 


C Mix 
0 


-Co 

Co -Co 


^ JftNc x 4 Nc 


\ c 0 -c 0 / 

Furthermore, again we assume to be in the reference situation described above, which implies that = {2} 
and = {1} for the domain Hi, i.e., the second coordinate of Ti is frozen. Let us consider a function 

M,2 1 J 

i/ 1 ) e V^ 1 (fti), which is written as 

v^\x) = ^ a^xp e (x), x € 

e£E 2 

Observing that 

v |r! (* c ) = (^( a oo^ + a oi^) tj>o(®i) + i a io a n i(^i))j 

we enforce = 0 if and only if the relations + c*oV ~ ® an< ^ + a i? = 0 are satisfied. In other 
words, t/ 1 ) has to be written as 

f (1) (z) = V2 (a'V^,o(ii) + (* 2 ), 

where r x is the univariate wavelet vanishing at 0, defined in (2.9). Since i/ 1 ) has already been set to 0 at 
C, we now enforce the linear combination 

~( a m + Q oV) + (“io* + «u) = 

The matrix containing the matching and boundary conditions along the sides takes the form 

e r N cxANc . 


Since there is no scaling function associated to C, we end up with the system 

Mix 

X> Mix a = 0, with Z> Mix = ^ ix eJR 2NcxANc . 

C\ 

Again, t> Mix is easily shown to have full rank, so that dim W^(il) = 2N C \ this is precisely the number of 
grid points surrounding C, to which the wavelets are associated (see Figure 5.3, center). 

Pure Dirichlet case. In case d), the matching and boundary conditions at C obviously coincide with 
those of cases b) and c). As far as the conditions at the sides are concerned, working out as before we end 
up with the matrix 
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and the whole system takes the form 


Z> Dir a = 0, with X> Dir = ° Djr g j*(2Jv c +i)x4Jv c _ 

Once again, D Dlr has full rank, so that dim Wf(Sl) = 2 N c - 1 (see Figure 5.3, right). 

5.2.1. Wavelets with minimal support. Let us first remark that, as for the interior cross point case, 
a wavelet supported in only one subdomain ft*, i = 2, . . . , N c - 1 has to satisfy (5.7). In the pure Neumann 
case, the extra orthogonality condition to <fj t c makes the existence of such a wavelet impossible. Conversely, 
in the mixed Dirichlet/Neumann- and pure Dirichlet cases, its existence is indeed possible. 

Pure Neumann case. Let us try to construct wavelets that are supported only in fii, 1 y V( „ , respectively. 
The corresponding conditions read 

(5.12) c 0 - a (1 ) = c' • a (1 ) =0, co • = c" ■ a^ c) = 0, 

which leads to « fl) = (a, b , -a, - b ) and a (/Vc) = (a, -a, b, -b). Adding the orthogonality condition yields, 
as before, a = 0. Choosing b= 1, we are led to define the following functions: 

t/,° c (x) := J^M*) “ 

L 0, elsewhere, 

i> 2N c c -\x) := { _ ^n{x), X e n Nc , 

1 0, elsewhere; 

in addition, let us define f° r ^ = — 1, i e {1, . . . , Nc~ 1} as in (5.9), ^ c for l = 2 i, i £ {1,..., Nc — 1} 

as in (5.10) and = ^j t c as * n (^.11). These are 2 Nc + 1 functions (the same number as the dimension 
of wj^fi)); two of them are supported in exactly one subdomain, one is supported in all the subdomains 
matching at C, while the remaining ones are supported across two consecutive subdomains. 

Let us now show that these functions are linearly independent. Indeed, let y’> := Y^m-o A Tn i'.’ T " r be such 
that t\) = 0. Now, ip(C) = 0 implies A 2 j v c = 0. On the other hand, V|r, = 0 implies A : = 0, and ib\r N = 0 
implies A 2 w c -i = 0. Thus, we are left with the same situation as in the interior cross point case. 

Pure Dirichlet case. As mentioned above, since in this case we do not have to enforce the orthogonality 
condition to <^,c, one can construct wavelets that are supported in only one subdomain. These are defined 
as follows 

(5.13) ^,2.-1^ ;= | V’oo(i) - V’oi(i) - ^io(i) +^n(x), xeSli, 

1- 0) elsewhere, 

for i g {1, . . . , N c }; they are precisely the functions introduced in (5.8) with a = 1. In addition, one has 
the wavelets , i € { 1 , . . . , N c - 1} defined in (5.10), which are supported in ff, Uf2 i+1 . No function 
supported in all the subdomains surrounding C is needed. 

We have defined a system of 2N C - 1 functions, which is precisely the dimension of The linear 

independence of these functions is readily seen. Indeed, since the linear combination ib := y^ 2 ^c ?— 1 \ im — 

* * Z_^J7l = l m ■ J,C 

0 must vanish on all sides, we obtain that the coefficients corresponding to the even indices arc zero. Next, 
considering all the subdomains, we get that the coefficients corresponding to the odd indices also vanish. 
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Fig. 5.4. Wavelets t jx* c and ip 3 c around a Dirichlet boundary cross point (Note again that only a portion of each 
subdomain around C is shown.) 


Mixed Dirichlet/ Neumann case. Again, no orthogonality condition to (p 3 ,c is needed. So, we use the 
system of functions constructed for the pure Dirichlet case, to which we add one function. Precisely, if 
Ih C r Neu ? we add the function 


/ 0 , X _ / ^io(x) -^oo(i), xetti, 

J ' C \ 0, elsewhere. 

If Tn c C ^ Neu j we add the function 

,2 N C ( \ f ^01 (X) - ^00(x), X € Qn C i 
J ' C \ 0, elsewhere. 

These functions are continuous and supported in fij, resp. (note that they correspond to the solutions 

cx( l ) = (-1,0, 1,0), q (Nc) = (-1,1, 0,0), resp., of (5.12)). As a whole, we have 2N C functions, so the 
dimension of Wf{$l) is matched. Finally, the linear independence is obvious. Indeed, assuming for instance 
I\ C r Neu , the Unear combination ip := A m ip™ c = 0 must vanisl1 on r x ; this implies A 0 = 0 and 

reduces the problem to the pure Dirichlet case. 


Example (continued). We consider an L -shaped domain made up by 3 square subdomains meeting at 
C . We enforce homogeneous Dirichlet conditions on the whole boundary. The two different types of wavelets 
produced by our construction are displayed in Figure 5.4 * ^ 

5.3. Matched wavelets across a side. Let us consider two subdomains and having a common 
side <7 := Pi Q_ . Moreover, let us denote by A and B the two endpoints of <r. As in the cross point case, 
we may reduce ourselves to a reference situation. Instead of thinking each subdomain as the image of the 
reference domain ft, here it is natural to think (Q+, resp.) as the image of the domain /_ x /+ (/+ x /+, 
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resp.), with 


(5.14) A = F_(0,0) = F + (0,0), B = F_(0, 1) = F+(0, 1). 

Since the reasoning for the reduction to the reference situation is analogous to the cross point case, we drop 
these arguments here. Now, (5.14) implies that the set of grid points on a is given by 


H a i — {h : h — F+(0, ^2 ) i ^2 £ Aj U Vj} 
— [h : h = F-( 0 , A2), /12 G Aj U Vj }. 


The grid points (Oj^j.i) and *+( 0 ,*'* Mj) are already associated to wavelets since these points 

correspond to the cross points A and B , respectively. Hence, we are left with the points 


K n< := n a \ {F + (0,v jA ),F + ( 0,v jlUi )}. 


Let us first consider the case h € Ji „ with /i 2 € A j. Consequently, there exists a scaling function 

{ fj, i(x_ + 1), if x = F_(x_,x 2 ) € fi_, 

0,o(x+), if 2 = F+(x + ,x 2 ) <E n+, 

0, elsewhere, 

associated to h. The basis functions of the local spaces Vh_ l (£l±) are then given by 


V*e (*) : =Sj y h 2 ( X 2) 


+ 1)» e — 0, 
(*-+!>, e-1, 


f Cj,o(i- 

\ ( 


£i,o(*+), e = 0, 
,(x+), e = 1. 


This shows that the matching at h is equivalent to the matching at a univariate interface point. Considering 
the wavelet functions ^ defined in (4.7) for l = — , +, we end up with the two wavelets 

^ ±(x) = 

1 0, elsewhere, 

which will be associated to the grid points /i_ := F-{uj ^ Afp ^ 2 ) and •= F+(z^i, / 12 ), respectively. 

Finally, we have to enforce the matching conditions along points h 6 where fi 2 6 Vj. In this case 

there is no scahng function associated to h. Again we are reduced to the univariate interface point case, but 
now in the situation considered in Subsection 4.2. Considering the basis functions tr- defined in (4.9) for 
l = — , 0, -f , we end up with three wavelets 

1 0, elsewhere, 


which will be associated to the grid points /i_ := F-(vj ho =: F-(l, fi 2 ) = F+ (0, / 12 ) and /i + := 

F+(v jy i,h 2 ). 


6. Tri variate matched wavelets. In this section, the exposition will be deliberately less detailed than 
in the two previous sections, as the three dimensional construction follows the same spirit presented before. 
We shall be mainly concerned with the matching around a cross point. The matching around an edge or 
across a face will be easily reduced to lower dimensional situations. 
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6.1. Matched wavelets around a cross point. Let us first assume that the cross point C belongs to 
ft. Let Nc = Nc t 3 denote the number of subdomains meeting at C, and let these subdomains be (re-)labeled 
by f2i, . . . ,Qn C ' It is not restrictive to assume that for all i e {1, ... , Nc }, one has C = F^(0, 0,0). 

It will be useful to express the number Nc , 2 of faces and the number Nc, 1 of edges meeting at C as a 
function of Nc ■ This can be accomplished as follows. Consider the tetrahedron in Cl 

T := {(xi,X 2 yX 3 ) : X\ > 0, x 2 > 0, x 3 >0, xi 4- x 2 + x 3 < 1} 

and, for each i (E {1, . . . , Nc}> set T x := F l {T). Then, Pc := Ui is a (distorted) polyhedron in JR 3 , with 
the following property: each face (or edge or vertex, resp.) of Pc is in one-to-one correspondence with a 
subdomain (or face or edge, resp.) meeting at C. Thus, by Euler's polyhedron Theorem, we get 

Nc , 3 ~ Nc , 2 + Nc, 1 = 2. 

On the other hand, since each face of Pc is a (distorted) triangle and each edge of Pc is shared by exactly 
two faces, one has 3iVc,3 = 2 JVc, 2 - It follows that 

(6-1) N c ,2 = § N c , N c , 1 = IJVc + 2. 

The grid points h eHj, to which we are going to associate the matched wavelets, have the form 

(6.2) hci,e = F,( C ei ,Ce 2 ,Ce 3 ), * € {1, . . . , N c }, e G E 3 \ {(0, 0, 0)}, 

with 

_ f 0, if ei = 0, 

Ce ‘ :_ W 1} if ei = 1. 

Note that there is exactly one of such points which lies inside each subdomain, each face and each edge 
meeting at C. Thus, the total number of such points is 

Nc , 3 4- Nc , 2 4- Nc , i = 3Nc 4* 2. 

This is precisely the number of wavelets to be constructed around C. Indeed, define for each e E E 3 

V’e’ ) ( x ) = i>e(x) = 0 ei (xi )6 e2 (x 2 )0 e3 (x 3 ), 

where 

_ f £j,o, if ej = 0, 

l Vj.Vj i i if e; = 1. 

Let us consider the spaces 

fj+i(Hi) == spanlV'e^ : e G £ 3 } = {v (<) = ^ : a (0 := (a^egE* G J? 3 } 

and let us introduce the column vector a (a^)^ € iR 8Nc (8 being the cardinality of E 3 ). In order 
to characterize the local space 

v£i(fl) := { V G C°(fi) : V|n. € if i G {1,...,JV C }, = 0 elsewhere}, 
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we enforce continuity firstly at C, secondly at the edges and finally at the faces meeting at C. This gives 
Nc — 1 conditions at C, which can be written as 


CqOL = 0, 


where 



-Co 


Co 


-Co 


G 1 R( n c-1)x8Nc 


with Co = (1, 1, 1, 1, 1, 1, 1, 1)- Next, we have N ed - 1 conditions at each edge, where N ed is the number of 
subdomains meeting at the edge; all these conditions can be written as C x a = 0 , for a suitable matrix C i 
whose structure depends on the topology of the subdomains. Finally, we have 1 condition at each face; they 
can be represented as C 2 ot — 0. 

In addition, we want to build functions in Uf(fi) := {v € Vf +l (Sl) : (v,^c)n = 0}; this adds the 
condition Ba = 0 , with B := (6, . . . , b) and 6 = (1, 0, 0, 0, 0, 0, 0, 0). Summarizing, we have 


T> ol = 0 


with 


T> = 


Co 
C i 

c 2 

B 


All the conditions that we have enforced are linearly independent, as shown in [4] (see Section 5.2); thus, 
their number is 


(Nc - 1) + J2 ( N ed - 1) + N c ,2 + l = N c + 3N C - N C ,1 + Nca = 5AT C - 2. 

edges 

Indeed, since each subdomain contains 3 edges meeting in C, one has 

( N ed - 1) = 3 Nc - Nc, 1 - 

edges 

We conclude that dim Wf(Q) = 3/V c +2, as desired. After the dual construction is made, biorthogonalization 
is accomplished as described in the previous sections; we omit the details. In conclusion, we end up with 

3Nc + 2 primaI and dual wavelets and c , l = 1, . . . , 3N C + 2; they are associated to the grid points 
(6.2), which from now on will be indicated by hc,i- 

Wavelets with localized support. Let us first notice that, as in lower dimension, no function in W c (fl) 
exists, which is supported in only one subdomain £li . On the contrary, it is easily seen that for each couple of 
continguous subdomains, two linearly independent functions in Wf(fi) can be built, which vanish identically 
outside the two subdomains. 

Thus, to each point hc,i which lies inside the common face of two subdomains, we associate one of such 
wavelet. To all but one points h c ,i lying inside the subdomains, we associate other such wavelets, choosing 
them to be linearly independent from the previous ones. To each point h c ,i lying inside an edge, we associate 
a wavelet supported in the closure of the union of all subdomains sharing the edge: it has the local structure 
of a tensor product of a two-dimensional scaling function (p j:C times the wavelet r in the direction of the 
edge. Finally, to the remaining point hc,i interior to a subdomain, we associate the global wavelet, which is 
the three-dimensional analog of (5.11). 
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Tensor products of matched functions. According to (6.1), the number of subdomains meeting at C 
is even. We consider here the particular situation in which these subdomains can be grouped in two sets 
of equal cardinality N c / 2, all the subdomains of each set sharing a common edge stemming from C. For 
instance, this is the relevant case of 8 subdomains meeting at C, and representing the images of unit cubes 
lying in the 8 octants of 1R 3 . 

By possibly introducing additional parametric mappings, we may reduce ourselves to the situation in 
which N c /2 subdomains lie in the upper half space x 3 > 0 while the remaining ones lie in the lower half 
plane. Each upper subdomain f li can be written as Q t = il[ x (0,1), where is a 2D subdomain in the 
plane x 3 = 0; the companion lower subdomain is D' x (-1,0). Thus, we are led to consider the case of N c /2 
subdomains in the plane meeting at C' = (0,0). Let tp 1 / be the bivariate scaling function associated to C ' , 
and let ( l = l,..., Nc/2) be any system of bivariate wavelets around C, built as in Subsection 5.1. 

Moreover, let (^ = — >”b) be the univariate wavelets defined in Subsection 4.1, and let dj (1 ,0, ) 

the univariate matched functions defined in Subsection 4.2. Then, a system of wavelets around C can be 
defined as follows: 

The functions , with l = l,...,iV c /2, are associated to the grid points having a strictly 

positive £3 -component and not lying on the £3 -axis; the functions ip* 1 an d are associated 

to the analogous grid points having zero or negative £3 -component. Finally, the functions p! 1 <8> i/jj are 
associated to the remaining grid points on the £3 -axis (obviously, these functions are extended by zero 
outside the union of the subdomains). 

In the case of 8 subdomains meeting at C , the bivariate scaling and wavelet functions may be chosen 
to be themselves tensor products of univariate matched functions, so one can obtain a fully tensor ized local 
wavelet basis around C. 

Cross points lying on the boundary . Let us now assume that the cross point C belongs to We follow 
the same notation as before. The grid points h € Hj around C to which wavelets will be associated are 
again of the form (6.2), but now the points lying on a face or an edge contained in f D j r are missing. Let us 
denote by resp.) the number of faces (edges, resp.) containing C and contained in f Dir . It is 

easily seen that the number lc of grid points hc y i we are interested in is 

lc = N C; . 3 + (Nc .2 - NS$) + (Nc.i - NS:i). 

Let us count the number of conditions that define the space Wj’ (fl). We have Nc,3 ~ 1 matching conditions 
at C, plus one vanishing condition if C £ foir or one orthogonality condition if C & foir- Next, we have 
N e d ~ 1 mat chin g conditions at each edge, plus one vanishing condition at each edge contained in C € Toir- 
Finally, let N j indicate the number of subdomains sharing the face / (this is 2 if the face is not contained in 
dn, 1 if it is); then, we have N f - 1 matching conditions at each face, plus one vanishing condition at each 
face contained in C G foir- Observing that 

^N ed = N f = 3Nc ’ 3 ' 

edges faces 

the total number tc of conditions which define (ffc) is 

tc = 7 Nc, 3 ~ (Nc , 2 - A£* 2 r ) - (No, 1 - ^\ r ). 

Since these conditions are linearly independent (see again [4]), we obtain, as desired, 

dim Wf(tl) = SN Cy 3 - tc = lc * 
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Fig. 6.1. Example of 5 trivariate subdomains meeting at a cross point belonging to the Dirichlet part of the boundary (left). 
The pictures in the right column show the wavelet grid points having positive, zero and negative z-component, respectively ( from 
top to bottom). 


As in the interior cross point case, wavelets can be constructed with localized support. Actually, if 
C £ r^jr, wavelets exist which are supported within one subdomain; they will be associated to the point 
he, i lying inside the corresponding subdomain. 

Example . Let us consider the domain represented in Figure 6.1 which is divided into 5 subdomains. 
Let us assume that the three boundary faces meeting at C are contained in Toir- We want to exhibit one 
particular choice of wavelets around C , which can easily be constructed. To this end , let us assume that 
C = (0, 0, 0) and let us divide the wavelet grid points around C into three sets , corresponding to their third 
coordinate z being negative, zero or positive as indicated in Figure 6.1. 

For defining the wavelets associated to the grid points having z < 0, consider the bivariate scaling function 
Tj and the bivariate wavelets ’ (l = 1, . . . , 6) associated to a 2D interior cross point C f common to three 
subdomains. In addition, let df be the matched univariate function defined in (4.9) and having support in 
(—1,0). Then , we associate to these grid points the 7 wavelets 

< Pj 1 ® (to the grid point on the z axis ) 

i {l = 1> • • • , 6) (to the grid points around the axis). 

It remains to define wavelets associated to those wavelet grid points having zero or positive third com- 
ponent. To this end , let ipj (l = 1, . . . , 3) be the system of wavelets associated to a 2D Dirichlet boundary 
cross point common to two subdomains. Then , the functions (l = 1,2,3) will be associated to the 

three wavelet grid points with z — 0 and the analogous three grid points in the upper half space are identified 
with the wavelets ipj 1 ’ 1 0 (l = 1,2,3). 

Note that no extra orthogonality to dual scaling functions has to be enforced , since no scaling function 
is associated to the 3D Dirichlet cross point C. □ 

6.2. Matched wavelets around an edge or a face. Let <r = ed be an edge, at which N e( i subdomains 
meet. We can reduce ourselves to the situation in which the subdomains are (re-)labeled by fli, . . . , Sl Ned 
and for each i e {1, . . . , N ed } we have 


6d={F i (0,0,C) : 0 < C < 1 }■ 
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It is enough to consider points h e d G £j+i which are internal to ed. Precisely, if h e d = ^(0, 0, ^ 3 ) with 
h 3 G A* nt , then we build wavelets which, in local coordinates, can be written as ^ ed < 8 > where ^ ed 

are matched bivariate wavelets as defined in subsections 5.1 or 5.2. They will be associated to the points 
h G Hj having the form h = F t (( 1 , £ 2 , h 3 ), with Ci, C 2 € { 0 , x}. 

On the other hand, if h e d = Fi(0,0,fi 3 ) with G V* n *, then the wavelets will be locally represented as 
'fij ed® 7 ! where now ed are matched bivariate wavelets, defined as in the previously quoted Subsections, 
but without enforcing the biorthogonality condition 23a = 0 . The association to the grid points surrounding 
h e d is done as above (note that now h e d G Wj). 

At last, let a — f be a face common to two subdomains and The reference situation is such 

that 


/ = {F + (0,<2,C3) : 0<C2,C3<1}. 


Let hj G /Cj+i. If hf = F+(0, & 2 i ^ 3 ) with h 2 , h 3 G Ay 1 *, then we build wavelets having the local represen- 
tation jL 1 t (g) £ t 0 f r , where f are matched univariate wavelets as defined in Subsection 4.1. On the 
other hand, if hj = F+(0, h 2 , h 3 ) with h 2 , fi 3 G A*- nt U VJ n * and at least one coordinate in V* nt , then the 
local representation of the wavelets will be one of the following ones: 






1 .h 


j,/l3 ’ 


PM 1 




where now tr ^ are matched univariate functions built in Subsection 4.2. The association of these wavelets 
to the grid points surrounding hf is straightforward. 


Appendix A. Mask coefficients. In this appendix we provide all the mask (filter) coefficients of the 
univariate scaling functions and wavelets for our B-spline example L — 2 and L = 4. These are the data 
that are needed to reproduce the figures in this paper and to use these functions, e.g., as trial functions for 
numerically solving differential and integral equations. Starting from the masks of the univariate functions, 
those for multivariate and matched functions can easily be obtained by appropriate tensorization, using the 
matching coefficients given in the paper at any interface among subdomains. 

Refinement coefficients . Let us start with the refinement coefficients on the real line corresponding to 
the equations 


= 2 1 / 2 ^ afc <t>{ 2z - /c), 

kez 


<j>{x) = 2 1/2 ^ dk 0(2x - /c), 
kez 


where the coefficients here are given by ([8]): 


a_i = 5 . OOOOOOOOOOOOe - 01, a 0 = 


l.OOOOOOOOOOOOe + 00, ai = 5. OOOOOOOOOOOOe - 01, 


for the primal scaling function, and by 


a_ 4 = 4.687 500000000e - 02, 

a_i = 5.937500000000e - 01, 

a 2 = —2.5 OOOOOOOOOOOe - 01, 


a -3 = — 9.375000000000e - 02, 

a 0 = 1.406250000000e + 00, 

a 3 = — 9 . 375000000000e - 02, 


a _2 = — 2 . SOOOOOOOOOOOe - 01, 

ai = 5.937500000000e - 01, 

a 4 = 4 . 687 500000000e - 02, 


for the dual one. 

Next, we give the entries of the refinement matrices Mj and Mj for the whole bases Ej, Ej, respectively, 
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as defined in (2.3-a). These matrices have the following structure (taken from [11]): 




Here, the central blocks A 3 and Aj take the form 

= 3 < m < 2 i+l — 3, 2<fc<2>-2, 

(A,) =~y=a m ^2k, 4 < m < 2 0+x -4, 4 < k < 2° - 4, 

with the previously given coefficients a k and a k . The upper left blocks are here given by 

" 7.07106781 1865e - 01 ' 

4.890821903207e - 01 

— 3. 794806392368e - 01 

— 6.010407640086e - 02 

M l = 

1.084230397819e - 01 
4.360491817317e - 02 

— 2.121320343560e - 02 
-1.060660171780c -02 . 


Ml= 


' 7.071067811865e — 01 
7.7339804 19228e - 01 
— 2.209708691208e - 01 
—3.3 1456303681 2e - 01 
1.657281518406e - 01 
0 . 000 OOOOOOOOOe 4- 00 
O.OOOOOOOOOOOOe + 00 
0 . OOOOOOOOOOOOe + 00 
O.OOOOOOOOOOOOe 4 00 
O.OOOOOOOOOOOOe 4 00 
O.OOOOOOOOOOOOe 4 00 


O.OOOOOOOOOOOOe 4 00 
7.402524115547e -01 
8.175922157469e - 01 
4.529902816976e - 01 
— 1.933495104807e - 01 
— 6.629126073623e - 02 
3.314563036813e - 02 
O.OOOOOOOOOOOOe 4 00 
O.OOOOOOOOOOOOe 4 00 
O.OOOOOOOOOOOOe 4 00 
O.OOOOOOOOOOOOe 4 00 


O.OOOOOOOOOOOOe 4 00 
— 2.695844603274e - 01 
— 6. 187184335382e - 02 
3.933281470350e - 01 
1.007627163191e 4 00 
4. 19844651 3 295e - 01 
— 1.767766952966e - 01 
— 6.629126073624e - 02 
3. 31 45630368 12e - 02 
O.OOOOOOOOOOOOe 4 00 
O.OOOOOOOOOOOOe 4 00 


O.OOOOOOOOOOOOe 4 00 
5.082329989778e - 02 
4.419417382416e - 03 
— 5.966213466262e - 02 
-1. 80091 2583334e - 01 
4. 19844651 3 295e - 01 
9.9436891 10436e - 01 
4.1 984465 13295e - 01 
— 1 . 767766952966e - 01 
— 6.629126073624e - 02 
3. 3145630368 12e - 02 


The lower right blocks Mr and Mr arise from their upper left counterparts by reflecting rows and columns 
in the following way: 


{Mr)2J — m,2J — k — ^Tl — 0, . . . , 8, Jc 0, 


and similarly for the dual functions. 

Wavelet coefficients . The wavelets on the real line can be written as linear combination of the translates 
of 4> on level 1, i.e., 

i>{x) — 2 _1/2 ^ b k <j)(2x - k ), 4>{x) = 2 _1/2 b k 4>( 2x - k). 

kez 
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The corresponding masks for our example are given by 

fci = 1.060660171780c -01, 62 = 2.121320343559c - 01, 63 = -5.656854249493e - 01 , 

64 = — 1 .343 502884 2 54e + 00, 65 = 3.181980515339e + 00, & 6 = -1.343502884254c + 00, 

bj = — 5.656854249492e - 01, b 8 = 2.121320343560e -01, 69 = 1.060660171780e - 01, 

for the primal wavelets, and by 

64 = -1.104854345604c - 01 , 65 = 2.20970869 1208e - 01, & 6 = -1.104854345604c - 01, 

for the dual ones ([8]). Let us now consider the two-scale matrices Gj and Gj for the whole wavelet basis, 
which give the transformation 

^3 ~ 1, T j = GjZ*j+\. 

Again, Gj and Gj have the same block structure as M 3 , M 3 , respectively, possibly with a different size of the 
blocks. The inner blocks correspond to those of Mj, Mj, by replacing a k , d k by b k) b k , respectively. Here, 
the upper left block of the refinement matrix for the wavelets read: 

7.07106781 1866e - 01 O.OOOOOOOOOOOOe + 00 ‘ 

— 4 . 89082 1 903 207e - 01 1.237436867076e + 00 

3.794806392368c - 01 -4.313351365238e + 00 

6.010407640085e - 02 4.985102807365e + 00 

— 1 .084230397819e - 01 -1.555634918610e + 00 

— 4.36049181 731 7e - 02 -6.717514421272e - 01 

2.121 320343560e - 02 2.121320343560e - 01 

1.0606601 71 780e - 02 1.060660171780e - 01 

7.07106781 1865e - 01 O.OOOOOOOOOOOOe + 00 

— 7.7339804 19228e - 01 - 1.4731 39 127472e - 02 

G l = 2.209708691208e - 01 -5.892556509888e - 02 

3.314563036812e - 01 1.325825214725e - 01 

— 1.65728 15 18406e - 01 -6.629126073624e - 02 

The lower blocks again arise by reflecting the upper blocks. 

Homogeneous boundary conditions. Let us now consider scaling functions and wavelets having homoge- 
neous boundary conditions. For the scaling functions, one only has to eliminate the first and last row in M 3 
and My, respectively. 

For the corresponding wavelets, we have to modify the first and last wavelet as defined in (2.9). This 
means changing the first and last column of Gi and (5^, respectively, by using the following coefficients 
instead: 


60 = 
63 = 
be = 

O.OOOOOOOOOOOOe + 00, 
8 . 500000000000e - 02 , 
3 . OOOOOOOOOOOOe - 02, 

bi = 
64 = 
b7 = 

— 6. 916666666667e - 01, 

— 1 .533333333333e - 01, 
1 . 5 OOOOOOOOOOOe - 02 , 

<§r ,? 

1! It 

5 . 366666666667 e - 01 , 
-6. 166666666667c - 02, 

for the primal wavelets, and 

bo = O.OOOOOOOOOOOOe + 00, 
53 = 4 . 68 7 500000000e - 01 , 

for the dual ones. 

bi = 
64 = 

— 1.093750000000e + 00, 
— 2.343 750000000e - 01, 

b 2 = 

3. 125000000000e - 01 , 


The coefficients in the refinement matrices for a whole variety of choices of the parameters L and L (for 
the construction in [11] based on [8]) can be obtained in MATLAB-format from the homepage of the third 
author http://www.igpm.rwth-aachen.de/~urban under the topic Software. The software described in 
[3] (written in C++) that has been used to produce the pictures in this paper can also be obtained via this 
website. 

Acknowledgement. The authors feel very grateful to Titus Barsch for his assistance in programming 
and producing the pictures. 
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